Color Coulomb Potential in Yang-Mills Theory from Hamiltonian Flows 
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A profound understanding of the confinement mech- 
anism in QCD still represents, after almost 40 years of 
intense research, one of the most important challenges in 
modern theoretical particle physics. We here report on 
a new approach to the subject which uses the recently 
developed technique of Hamiltonian flows The gen- 
eral setup is the Hamiltonian formulation of Yang-Mills 
theory in the Coulomb gauge 0. Important progress 
has been made over the last decade in this formulation, 
mainly via the variational principle [3l4l2j. The hori- 
zon condition is implemented in accord with the Gribov- 
Zwanziger confinement scenario (l3l . IbH ] , and scaling be- 
havior of the equal-time two-point correlation functions 
results together with an infrared fixed point of an appro- 
priately defined running coupling constant (loj . 

The potential between static color sources in the large- 
distance or infrared limit is usually approximated by the 
color Coulomb potential. The color Coulomb potential 
is given by the vacuum expectation value 



((-dD)-\-d 2 )(-dD)-i) 



(1) 



[see Eq. (U) below for our definition of the covariant 
derivative D] and is usually expressed as 



((-dDy^i-d^fi-d^d-dD)- 1 ) 



(2) 



with the so-called Coulomb form factor /, which satisfies 
a Dyson-Schwinger equation (DSE). In order to calcu- 
late the color Coulomb potential, in Ref. the Coulomb 
form factor was simply set equal to one, while in Rcfs. 
0, HI, [HI the DSE for this form factor was approximated 
by replacing in the loop integral the full ghost propaga- 
tor dD)^ 1 ^ with the bare one, which results in an in- 
frared finite Coulomb form factor. In this way, a strictly 
linear growth of the color Coulomb potential with the 



distance between the color sources (for sufficiently large 
distances) has been found in Ref. [1J|. It would now be 
natural to try to improve the approximation by using the 
full DSE for the Coulomb form factor. However, it turns 
out @, H] that the full DSE for the Coulomb form factor 
cannot be consistently solved together with the DSEs for 
the static (equal-time) gluon and ghost propagators with 
an infrared-divergent ghost form factor, i.e., implement- 
ing the horizon condition. In other words, a confining 
color Coulomb potential cannot be obtained within the 
present approximation if the full DSE for the Coulomb 
form factor is used. 

In the present letter, we will focus on the determina- 
tion of the color Coulomb potential with the help of a 
different functional technique, the Hamiltonian flows [l[ . 
Interestingly, a consistent solution which exhibits scaling 
behavior of the static propagators and the color Coulomb 
potential is readily found in this framework, without any 
additional approximation for the Coulomb form factor. 

The organization of this letter is as follows: we will 
start with a brief presentation of the Hamiltonian flow 
technique and summarize the results of Ref. [l[ . We then 
derive the flow equation for the color Coulomb potential 
and, finally, present and discuss its solution. 

The construction of the Hamiltonian flows starts from 
the generating functional for Green's functions at equal 
times, 



Z[J] = / VAdet(-dD)\ip[A] 



2 e J - A 



where D denotes the covariant derivative 
Df = S ab di - gf abc A1 



(3) 



(4) 



and the dot in J • A stands for the contraction of color 
and spatial indices and the integral over position or mo- 
mentum, 



J- A 



d 3 xJf(x)^(x) 
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d 3 p 
(2tt)3 



J?(-p)A 4 (p). 



(5) 
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In Eq. (|3|), ip[A] represents the vacuum wave functional, 
det(— 3D) is the Faddeev-Popov determinant in Coulomb 
gauge 0, and the functional integral is over the trans- 
verse spatial gauge fields that fulfill the Coulomb gauge 
condition diAf = (we denote the contravariant spatial 
indices as su&indices). 

We now write the Faddeev-Popov determinant in the 
well-known way as a functional integral over ghost fields 
and introduce additional source fields for the latter. Fur- 
thermore, we introduce an infrared regulator term ASk 
with infrared cutoff scale k. Then the fc-dependent gen- 
erating functional becomes 

Z k [J, a, a] = J VAVcVce- s ~^ Sk+J - A+ °- c+e -' J (6) 

with the "action" 

S = -ln|^[A]| 2 + J d 3 xc a (x)(-dD) ab c b (x). (7) 



The regulator term is 



AS k [A,c,c} = \ J J^La?{- P )R a 



;(p)aUp) 



-0^c a (-p)gR c Ap)c a (p), (8) 



with the regulator functions chosen in the present work 

as 

(k 2 v 2 \ 
R A>k (p)=2pexpl^-^\ , 



' k 2 p 2 

RcAp) = P 2 ex P ( - p 



(9) 



Here and in the following, we use the notation p = |p|. 
The change of Zk as defined in Eq. ([5]) under a change of 
k constitutes a (functional) renormalization group trans- 
formation. 

In the limit of large k, the functional integration in 
Eq. §Q is restricted to ultraviolet modes where asymp- 
totic freedom sets in, so that Zk can be evaluated in per- 
turbation theory. In the other extreme, for k — >• 0, the 
cutoff term ASk vanishes and Zk — » Z, the full generat- 
ing functional. Between these limits, the change of Zk 
with k obeys an exact intcgro-diffcrcntial equation, the 
functional renormalization group equation or flow equa- 
tion. 

In practice, it is more advantageous to consider the k- 
dependent effective action Tk, the generating functional 



of one-particle irreducible Green's functions, instead of 
Zk- We then approximate the corresponding exact flow 
equation for the variation of Tk with k, by approximating 
the form of Tk- Our basic assumption is that there exists 
a solution which shows scaling behavior in the infrared. 
We parameterize the static propagators as 



(27T) 6 



6 2 \nZ k 



5J?(-p)5J (<D 



J=a=S=0 



G A , fe (p)rt u (p)(2 7 r) 3 5(p-q). 
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6 8 2 lnZ k 



Sa a (-p)S<r b (q) 
1 



J=a=S=0 
ab /r* \3 



fJ 



G c ,k(p)S ab (27rrS(p-q), (10) 



with the functions 

GaMp) 

G c ,k(p) 



1 



2oj k (p) + R A ,k(p) ' 
1 

p 2 /d k {p) + RcAp) 



(11) 



In the first of Eqs. (fTO]). ijj(p) denotes the transverse 
projector or spatially transverse Kronecker delta. Our 
assumption of scaling then refers to a power behavior for 
small momenta of the functions LOk=o{p) and dk=o(p), 



u) (p -> 0) oc p 



d (p -> 0) cx p 



(12) 



The scaling assumption has a number of important 
consequences, just as its counterpart in Landau gauge. 
As shown in Refs. [I, @, [13, [3 [II , the exponents fulfill 



a = 2/3-1, a>-3/2. 



(13) 



In particular, in the infrared the diagrams with the great- 
est number of ghost loops dominate over the others (ghost 
dominance). Furthermore, the dressing of the ghost- 
gluon vertex is constant [IH, [lj|, usually referred to as 
the nonrenormalization of the vertex (see also Ref. [20|). 
so that we can calculate with a bare ghost-gluon vertex, 
possibly adjusting the value of the coupling constant. Fi- 
nally, we can neglect vertices with more than two gluon 
lines, as far as the dominant infrared behavior is con- 
cerned. 

Retaining only the contributions that are relevant to 
the infrared behavior, the flow equation for Tk induces 
the following intcgro-diffcrcntial equations for the static 
propagators: 



J 



a n c 

dk^^-- 



d 3 q 
(2tt) 3 



- dR^k- 

^rc,k — -q^ — "ci 



( 9 )G c , fe (|p + q\)q 2 (1 - (p • q) 2 ) , 
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Ok 



d- k 1 (p)=N c 



d 3 q 
(2^ 



..fc^r^Gx,* )(?)Gc,*(|p + q|) 



G, 



dk 
' dk 



G c , k (g)Gx, fc (|p + q|) 



(p + q) 



(l-(p-q) 2 ) 



(14) 
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FIG. 1: Flow equation for the gluon propagator. The spiral 
and dotted lines with the black circles denote the full gluon 
and ghost propagators at cutoff momentum k, respectively. 
An insertion of the fc-derivative dRk/dk of one of the regulator 
functions is represented by a square with a cross. In this and 
the following figure, the diagram on the left-hand side stands 
for the inverse propagator with the regulator function omitted. 

FIG. 2: Flow equation for the ghost propagator. 

where we have neglected in each equation the contribu- 
tion of a tadpole term, in order to obtain a closed system 
of equations. The Eqs. (|14p are represented diagrammat- 
ically in Figs. Q] and [2] For details of the derivation of 
these equations, we refer the reader to Ref. [l[. 

The solution of Eqs. (|14p is determined by imposing 
initial conditions on LJk(p) an d dk{p) at k = A, where A 
is an arbitrary large momentum scale. Here we adjust 
these initial conditions to satisfy the normalization con- 
ditions that u}q(p) oc p for large momenta p < A, and 
that rfo(p) follow a power law in the infrared (p — > 0), 
reflecting the properties of asymptotic freedom and in- 
frared scaling, respectively. It is possible to fulfill these 
normalization conditions by adjusting the constants a in 
w a(p) = a + p and d^{p) = g?a, technically upon rewrit- 
ing the differential equations ([M]) as integral equations 
over k. Again, we refer the reader to Ref. [l[ for details 
of the numerical procedure and here only cite the result, 
represented in Fig. [3l Note that for technical reasons, 
the equations can only be integrated down to a minimal 
value £; m i n of k, and we show the dependence on the value 
of fc m in in the figure. 

In the infrared, for p > fc m j n , we find the power behav- 
ior of Eq. (|12p. with the numerical values for the expo- 
nents 

a = 0.28, ,3 = 0.64, (15) 

which conform to the sum rule (|13[) . a = 2/? — 1. Note 
that a power behavior in the infrared (without specifying 
the exponent) was implemented as normalization condi- 
tion for the ghost dressing function, but not for the gluon 
propagator. We also find an infrared fixed point for the 



(running) effective fine structure constant, 

47T UJ (p) 

We can compare these results with the solutions of 
DSEs obtained in the variational approach in Rcfs. 
0, H, El- In this approach, two different solutions arc 
found with ft = 0.80 @, and ft = 1.00 [n| and the 
respective values of a according to the sum rule (|13l) . 
The same values of the exponents are also obtained an- 
alytically in an approximation that is only adequate in 
the infrared region ((J [Io[ ■ It is not unexpected that the 
functional renormalization group yields smaller values for 
the critical exponents than the DSEs in corresponding 
approximations. In similar calculations in Landau gauge 
Yang-Mills theory it has been found that the precise val- 
ues of the exponents vary slightly with the choice of the 
regulator functions. The upper limit of the interval over 
which the values of the exponents vary, is given by the 
result of the corresponding DSEs and is obtained from 
the functional renormalization group for so-called "opti- 
mized" regulator functions (2lj . 

In fact, we can simulate the use of optimized regu- 
lator functions by simply replacing the functions iOk(p) 
and dk(p) on the right-hand sides of Eqs. (fTl)) by u>o(p) 
and do(p). We have argued in Ref. pj that this replace- 
ment furthermore effectively takes into account part of 
the tadpole terms that we have omitted so far. With 
u>o(p) and do(p) on the right-hand sides of Eqs. ([T4"|) . we 
can perform the integration over k analytically and end 
up with equations very similar to those of the variational 
approach. We have shown in Ref. jlj that the numerical 
solution of these equations matches almost perfectly the 
one of Refs. @, H[, with (3 = 0.80. In this sense, our so- 
lution (fT5)) corresponds to the solution (3i = 0.80 of the 
variational (DS) approach. 

We now come to the calculation of the color Coulomb 
potential in the functional renormalization group ap- 
proach. For its definition, one considers the theory in 
the presence of an external static color charge density. 
Then the color Coulomb potential is the vacuum expec- 
tation value of the part of the Hamiltonian that depends 
on the external color charges, explicitly in momentum 
space (in an integral kernel notation) 

^ Q V-q) = ((p,a|(-5 J D)- 1 (-5 2 )(-^)- 1 |q,6)) 

= V c (p)6 ab (27r) 3 5(p- q ). (17) 

Introducing the composite operator 

K = J (0? c " d (- k ) fcV ( k )' 
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FIG. 3: Double-logarithmic plots of the inverse gluon propagator Wfc(p) (left) and the ghost dressing function dk(p) (right) 
resulting from a numerical solution of the flow equations (|14|) . for different minimal cutoffs k = fc m i n - The scale on the 
momentum axes is an estimate obtained by matching the position of the minimum of the inverse gluon propagator to the one 
in Ref. El. 



one can write 

F afc (p,-q) = ( C a (p)Xc b (-q)) GI . (19) 

The label GI on the vacuum expectation value stands 
for gluon-irreduciblc, meaning that one has to restrict 
the contributing diagrams to those where the operator K 
remains connected to the external points when all gluon 
lines are cut. 

The fc-dependent color Coulomb potential F k b (p, — q) 
is then naturally defined by including the cutoff term 
ASk in the functional integral representation of the vac- 
uum expectation value (fT9| as in Eq. ([6]), and a flow 
equation for F k can be derived in the standard way. 
For reasons of space, however, here we present a much 
quicker and equivalent derivation of the flow equation 
for F k which is based on the identity [22| 

^ [g(-dD)- 1 ] = (-dD)-\-d 2 ){-dD)-' (20) 

for the operators. We generalize this identity to 

|- [g{-dD + gR c>k )- 1 ] 

= (-8D + gR c , k )- l {-d 2 ){-dD + gRc.k)- 1 (21) 
for our present purposes, so that 

*f (p,-q) = [9(-9D + gRck)- 1 ] M» fc 

= V c , k (p)S ab (2n) 3 6(p-q). (22) 

Note that the rescaling of the ghost regulator function 
with a factor of g is essential to achieve a form equivalent 
to the definition of Fk described above. 



In order to put the identity (|2"Tj) to use inside the vac- 
uum expectation values we arc interested in, we define a 
^-derivative "at fixed integration measure" 

d 9 \ fm (°)k = J VAdet(-dD + gR c>k ) (i°i A ]) 

x \ip[A}\ 2 exp(-^A-R Ak -A\ , (23) 

for an arbitrary operator This definition immedi- 

ately implies that 

d g \ im G A , k =0. (24) 

For the application of the ^-derivative to the static ghost 
propagator, we use the identity 

(c a (p)c b (-q)) k = {(p,a\(-dD + gR c ^- 1 \ q ,b)) k . 

(25) 

With the help of definitions (JTOJl and Hfity, we then find 

d g \ im G C!k (p) = V c , k (p). (26) 
Introducing the Coulomb form factor f k (p) by 

V c . k (p) = \G c , k (p)p 2 f k ( P )G c , k (p) , (27) 

g 2 

we may rewrite the latter identity as 

d 9 \ im d-\p) = 1 d g \ !m G-l( P ) = ~h(p). (28) 

As a consequence of these relations, we can derive a 
flow equation for the Coulomb form factor by simply dif- 
ferentiating the flow equation (fT4| for d^ 1 with respect 
to g with the result 
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m h<p) = -N e 



d 3 q 
(2ir) 3 



J^^Cu, ) (q) G 2 c J\p + q|)(p + q) 2 / fe (|p + q| 



Ok 



dR 



+ 2 ( G c , k ^Gl k ) (?) q 2 f k (q)G A M\p + q|)- 



p + q) 5 



fl-(p-q) 2 ) , 



(29) 



where we have used the fact that d g \ { and the k- 
derivative commute. Since we have derived Eq. (|29[) from 
Eq. (|14[) . several approximations are implicit in Eq. (|29[) 
corresponding to those employed before in the derivation 
of the flow equation for d7 . 

The standard derivation of the flow equation that 
makes use of the composite operator K [defined in Eq. 
(|18[) ] and its equivalence with the argument presented 
above, as well as an algebraic construction that avoids 
reference to gluon-irreducible diagrams as in Eq. (fH?)) will 
be detailed in a future publication. 

Since the flow equations (fT4")) for Uk{p) and dk(p) do 
not involve fk(p), we can insert the solutions of the latter 
equations found in Ref. [l| in Eq. (f2T)|) and integrate this 
flow equation applying the same techniques used before, 
i.e., we convert Eq. (|29[) into an integral equation and 
solve it numerically by an iterative procedure. Equation 
(j29|) is linear and homogeneous in fk(p), and we decided 
to normalize fk{p) to one in the ultraviolet (below the 
initial scale A) by appropriately adjusting the initial con- 
dition f\(p) = /a- It is clear from perturbation theory 
that fk{p) should be constant in the ultraviolet (as long 
as k <C p) except for logarithmic corrections. 

Somewhat surprisingly, and contrary to the negative 
result of the search for a scaling solution in the variational 
approach complemented with DSEs , a solution of Eq. 
(f2T)| is readily found in the way described and is here 
represented in Fig. With the propagators obtained 
from the flow equations (|T4"|). we get a power behavior 

fk( P )^p-' 1 (30) 

in the infrared, with 7 = 0.57 for p > fc m ; n - Consequently, 
the color Coulomb potential behaves as 

V Ctk (p) cx p~ s , 5 = 2 + 2(3 + 7 (31) 

for p > fc m i n , see Eqs. (fTTj) and ([27|) . Making use of our 
result ([15]) for /?, we extrapolate V c ,k{p) to 

V c (p -> 0) oc p- 3 85 (32) 

for k = 0. We thus come quite close to a (p -4 )-behavior 
which would correspond to a potential that rises lin- 
early with distance (for sufficiently large distances) . Also 
note that our present result 0.64 for (3 is supposed to 



be smaller than the correct value as discussed below Eq. 
(|16l) , hence an improvement of the current approximation 
is expected to enhance the infrared exponent of V c (p). 

The result (|3"2"j) has been obtained with the propagators 
taken from the flow equations (TTJJ that do not include 
the tadpole diagrams. As argued before below Eq. (|16|) . 
we can easily take a part of the tadpole contributions 
into account in an effective way by replacing LUk{p) and 
dk(p) on the right-hand sides of Eq. (fT4"|) with u)o(p) and 
do(p). In Ref. [l| we have also computed the flow of 
this improved truncation leading to /3 = 0.80. Inserting 
the latter solution into the flow equation (|2U)l for the 
Coulomb form factor, we are led to Fig. @] (right). The 
resulting infrared potential reads 

V c (p-+0) ocp- 4 - 25 . (33) 

The exponents in Eqs. (|32p and Q33p provide us with an 
estimate for the systematic error of the present approxi- 
mation: 

V c (p -> 0) oc p- s with S e [3.85 , 4.25] , (34) 

including <5 = 4. 

In summary, the method of Hamiltonian flows allows 
for scaling solutions for the static two-point functions and 
the color Coulomb potential without additional approx- 
imations for the Coulomb form factor, contrary to vari- 
ational approaches. In particular, we find an infrared- 
divergent Coulomb form factor and an almost linearly 
confining potential V c (p) oc p~ s with S E [3.85, 4.25]. 
We expect that an improvement of the approximation 
employed narrows the above interval for 6 while still in- 
cluding 5 = 4. 
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